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POD-based Reduced Order Modeling 


Clarkson 


UNIVERSITY 


► Develop low-dimensional model that can represent the 
dynamics of a higher dimensional system 

► Uses: 

► Models for real-time control 

► Sub-model generation 

► Difficulties: 

► Round-off sensitive 

► Unstable / poorly conditioned 

► Units inconsistency for systems 
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Clarkson 


Particle Modeling 


UNIVERSITY 


► Incompressible, flow over a sphere with variable inlet velocity 

u = Uoc (1 + As\n(27rt/T)) 


► Axial Velocity 



► Centerline pressure 



► Force 
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Generate Modes 


Clarkson 


UNIVERSITY 


► 

► 


Find decomposition ^ - aj(t)4>j{x) to represent solution ty(x, t ) 


Choose optimal functions by maximizing mean square 
projection: 




((6 uf) 


► Some trouble: 



u r 

u z 


P 






Dimensionally inconsistent? 
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Possible Choices for u 


Clarkson 


UNIVERSITY 


► Depends on non-dimensionalization: 

U = ( U r (x , t), U z (x, t ), p(x, t)/(pUoo)) 

► Could be imaginary: 

u = (u r (x, t ), u z (x, t ), \/2(p(x, t))/p) 

► Also has arbitrary constant: 

u = (u r (x, t), u z (x, t), v / 2(p(x, t) + po)/p\ 
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Clarkson 


Eigenvalues of Steady Problem 


UNIVERSITY 



► 31 steady snapshots from Re = 0.1 to 100 

► Eigenvalues exponentially decay 

► Small number of modes can capture most of the energy 

► Last five modes of Lapack DGESVD are negative 
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Generation of Reduced Order Model 


Clarkson 


UNIVERSITY 


Lh 

o 


► Failure of Galerkin Projection 



T& + (V- V)\/ = -Vp + //V 2 v 

V • v 


dn = o 


► Modes are all incompressible 

► Continuity is always zero 

► Vp term integrated by parts is zero 

► Pressure modes may be undetermined 
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Clarkson 


SUPG Projection 


UNIVERSITY 




T 


c/fi + / (j) e, f • nrc/r 


+ 



dcj) de 
du Jr 


d(j) di 
dr] du 




► Streamwise-upwind-Petrov-Galerkin variational approach 

► Allows us to seek u n u z ,p with no pressure decoupling 

► Upwinds all the terms in the governing equations consistently 

► Results in a system of M ODE’s (Solved using DIRK & 
Newton-Rhapson) 
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Clarkson 


Drag Results - Steady 


UNIVERSITY 



► Accuracy of drag prediction versus Re for 5 and 10 modes 

► Significant improvement over DNS (10 vs. 21,000 degrees of 
freedom) 

► Close to empirical correlations in performance 

► All three sets of modes perform similarly 



NASA/CP— 2010-216112 504 


Clarkson 


Unsteady Data 


UNIVERSITY 


► Re = 0.1, 1, 10, 100 x St = 0.1, 0.5, 1, 2, 10 

► 15 time-steps/period, 20x15 snapshots 



► All eigenvalues after « 60 are probably round-off dominated. 

► Using Lapack DGESVD, all modes beyond 180 have negative 
eigenvalues 
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Unsteady Results 
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Clarkson 


UNIVERSITY 


Re = 0.1, St = 10 



Re = 100, St= 10 
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Clarkson 


Comparison of Modes 


UNIVERSITY 


Re = 1, St= 1 



► Re = 1 and St = 1. 

► SVD modes seem ok in spite of eigenvalue distribution. 

► Low Mach # modes give similar results with 20 modes, but 
for 60 modes Newton diverges. 
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Clarkson 


Conclusions 


UNIVERSITY 


► The devil is in the details!!! 

► Need stable numerical methods 

► Round off error can be considerable 

► Not convinced modes are correct for incompressible flow 

► Nonetheless, can derive compact and accurate reduced-order 
models. 

► Can be used to generate actuator models or full flow-field 
models 



